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Abstract 

We prove that when the bag surface is allowed to move radially, the equations 
of motion derived from the MIT bag Lagrangian with massless quarks and 
a spherical boundary admit only one solution, which corresponds to a bag 
expanding at the speed of light. This result implies that some new physics in- 
gredients, such as coupling to meson fields, are needed to make the dynamical 
bag a consistent model of hadrons. 
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The MIT bag model, in which hadrons are modelled by the states formed with free quarks 
confined inside an impenetrable bag, has been rather successful in reproducing static ground 
state properties of hadrons Because of its simplicity, especially its spherically symmet- 
ric version, the model has been used extensively in the discussion of various phenomena 
ranging from strange stars to ultra-relativistic heavy-ion collisions , even though these 
often involve situations of high density/temperature where the applicability of the model 
is doubtful. Many attempts have also been made to augment the basic MIT bag model 
with new ingredients, such as the partial restoration of chiral symmetry via meson coupling 
Q and inclusion of perturbative gluon exchanges among quarks However, in almost 
all discussion, the hadron bag is treated as a static boundary between the perturbative and 
nonperturbative vacuum, and excitations of the hadron are associated with the quark degree 
of freedom. The few notable exceptions [P-[10[|, which allowed for the possibility of a dy- 
namical bag boundary, focused mainly on reproducing the correct phenomenological parity 
order of the low-lying states of the nucleon, but several approximations and modifications 
to the theory had to be employed. Moreover the more fundamental question of whether it is 
consistent and feasible to use a dynamical bag to model hadrons was not addressed. That 
is the motivation of the present work. In this paper we discuss, analytically and without 
approximations, the consequences of allowing the bag boundary to move. 

We consider the basic MIT bag model with free fermions inside a spherically symmetric 
but non-static bag. We will show that the equations of motion require that the fermion 
field ip{t,r) vanishes at the bag boundary r = R{t). Hence our problem reduces to a quite 
general one: that of the Dirac equation in a spherical dynamical cavity with %p{t,R{t)) = 
0. Furthermore, we will show that for massless fermions, the only solution is that of a 
bag expanding with the speed of light. From this unexpected result, which evidently has 
no phenomenological application, we infer that some new physics ingredients, such as an 
interaction term with other fields, have to be introduced to make the more general dynamical 
bag model consistent and physical [|n|]. In the case of massive quarks we could not find a 
solution and we conjecture that in fact it does not exist. 
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The MIT bag Lagrangian density is written as 



where 9y{x) is 1 inside the bag and outside and 
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being the surface delta-function and the normal vector to the bag. From the Euler 
Lagrange equation of motion 

dC 
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we obtain 



i^'^d.ilj = 



inside the bag 



(4) 



on the surface 



(5) 



This last equation may be considered as the boundary condition for Eq. ^. Energy- 
momentum conservation implies a further constraint at the boundary p|,p!2|: 



lid' {^] 



r=R 
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We now look for a spherically symmetric solution of the above equations. In this case 
we have 9y = 9{R — r), = 6{R — r), and 



= {R, -r) . 



We first find an explicit expression for the boundary condition. Eq. H becomes 
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with 
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We can write the spinor as fT^JT^ 

fA ( 



(10) 



where y^f contains the spin and angular parts of the wave function. Observing that [a 
r)yji = —yfi' mil) we can write Eq. ^ as 



zRg{t,R)-f{t,R) 
Rf{t,R)-tg{t,R) 



g{t,R) , 
tfit,R) . 



(11) 



If -R = we have the famihar boundary condition for the static MIT bag, i.e. g{t,R) = 
—fit, R), which also corresponds to the one in the Bogolioubov model [|1^], and the analytical 
solutions are well known. However it is easy to verify that, if R ^ 0, then Eqs. |Tl| can be 
satisfied only if 



git,R) = f{t,R) = 0. 



(12) 



Notice from Eq. ^ that this implies B = 0, and in this way energy-momentum is conserved 
regardless of the motion of the bag. This also means that Eq. |^ does not provide information 
about R. 

Performing the change of variable y = rRQ/R{t) one could recast the problem into a 
static boundary one. In this framework the motion of the bag's surface is treated as a time 
dependent perturbation to the static Hamiltonian and one looks for the solutions by means 
of a time dependent expansion in terms of the static cavity eigenfunctions, exp{—iEnt}il'n{y) 



T5| . However in our case this approach cannot provide us the solution. In fact, writing the 



wave function tpf for the fixed boundary problem as 

oo 

^/ = E c„(t)e-^^"Vn(?/) 

n=0 



(13) 



we can in general work out c„(t), e.g. by perturbation theory, but it is well known that 
the static eigenfunctions, exp{—iEnt}ipn{y) , of the Dirac field inside a spherical cavity are 
non-zero at the boundary |]12[, and we have no way to impose the boundary condition Eq. |12 



on expression Eq. |T3|. In other words, although at the initial time we can always choose a 
suitable combination of ipniv) which is zero ai y = Rq, at subsequent times ipf(t,Ro), as 
expressed in Eq. |1^, will in general be different from zero. We hence need to proceed in a 
different way. 

Substituting expression Eq. ^ for ip in Eq. ^, with / = and /' = 1 in order to have 
spherical symmetry, we obtain two coupled equations: 



dt dr r 



■ df dg 



(14) 



(15) 



dt dr 

Integrating Eq. |l^ with respect to time and substituting the expression for f(t, r) in Eq. 



we have 



(16) 



.dg_ 

''dt 



dt 



d'^g 2dg_ 

Q^2 J, 



dA{r) 2 ^ , , 
dr r 



(17) 



where A{r) is a time independent function to be determined. Eq. [T^ is the spherical wave 
equation in integro-differential form, whose general solution can be written in the form: 



g{t, r) = - \G{t - r) - G{t + r)] + Ct + D 
r 



(18) 



with C and D two constants. In principle, a term of the form a/r is allowed in Eq. 18. 
However, it makes the wavefunction not normalizable and unphysical. We have also checked 
that this term does not affect at all our proof. We therefore set a = 0. 
Inserting the above expression in Eq. [T^ yields an equation for A{r): 



dA 2 ^ 

— = — A - iC 

dr r 



whose solution is 



(19) 



A{r) 



—Cr . 
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(20) 



Hence we can write an explicit expression for / as 

fit, r) = -{Git -r)+Git + r) + - [Qit - r) - Qit + r)]} - t^r , (21) 
r r 3 

with Q'iz) = Giz). Now using the boundary conditions Eq. |l^ we can derive the relations 
among Git — R), Git + R), and Equating git, R) and fit, R) to zero we have 

git, i?) = 4 [Git -R)-Git + R)]+Gt + D = 0, (22) 
R 

fit, R) = ^{Git -R) + Git + R) + hqit -R)- Qit + R)]} - i^R = . (23) 
R R o 

Taking the time derivative of both equations we obtain 

j^git, i?) = 1 {(1 - R)G'it -R)-il + R)G'it + R)} - ^ [Git - R) - Git + R)] + G = 

(24) 

and 

if it, R) = -|/(t, R) - j^qGRR - (1 - R)G'it - i?) - (1 + R)G'it + R)+ 



(25) 



Ji. 

i?2 



[Qit -R)- Qit + R)] + ji [{R - l) Git - i?) + (i? + l) Git + R)]} = 0. 



From the last two equations, and by means of Eqs. ^ and |23|, we can find 

(1 - R)G'it - i?) - (1 + R)G'it + R) + RiGt + D) + GR=0 (26) 

and 

- (1 - R)G'it - i?) - (1 + R)G'it + R) + Gt + D + GRR = . (27) 

These two equations contain all the information we need to solve the problem, i.e., to find 
i?(t) for any t > t^ and Giz) for any z > t^ + -R(to)- Iii fact, using the same argument as in 
|T^, as long as |-R| < 1, at each time t > to G'it — R) is known and we have two equations 
for the two unknowns Git + R) and -R(t). By summing and subtracting Eqs. ^ and we 
can decouple R and G'it + R) as follows 
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(l + r) [-2G'{t + R) + {Ct + D + CR)] = , (28) 

(l - r) [2G'{t ~R)-{Ct + D- CR)] = . (29) 
It is finally evident that Eq. implies 

R{t) = 1 , (30) 

and from Eq. ^ we have 

G' (t + Rit)) = ^[Gt + D + CRit)] . (31) 

Notice that -R = — 1 is not allowed because Eq. |29| would not be satisfied. Since ex- 
pressions Eqs. |18| and ^ represent the general solution of the problem, the above result 
excludes the possibility of any other solution. It is also important to note that, with = 1, 
t-R{t)=to-R{to). 

Now, knowing that R{t) = Rq + t — to and defining z = t + R{t), we have 

G'{z) = ^{Cz + D) z>to + Ro. (32) 

At this point, in order to have a clear understanding of the solution, we set C = D = 0. In 
fact, the two constants, with the boundary conditions given by Eqs. and are physically 
irrelevant. We obtain 

Giz) = Gizo) = G{to -Ro) z>to + Ro = zo, (33) 

Q{z)=G{zo){z-Zo)+Q{zo) z>to + Ro = zo. (34) 

Given ilj{to,r), which for example is different from zero at some r, then ip{t,r) will go to 
zero as t ^ to + Ro + r , and so we will have 

^p{t, r) = , r <t-to- Ro , r > t - to + Ro ■ (35) 

Therefore, after a time t = Ro + to, the solution represents an expanding spherical shell 
with internal radius Ri^ = t — to — Ro and external radius -Rout = t — to + Rq. The solution 
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for t < to can be found analogously, evolving backward in time. In this case Eqs. ^ and ^ 
imply R = —1 and G'{t — R) would be the unknown. Obviously the solution becomes 
singular as i?(t) = 0. 

This unexpected solution of the massless Dirac equation in a spherical dynamical cavity 
is evidently due to the boundary conditions Eq. However, as far as one considers the 
Dirac field completely confined in a cavity, any other boundary condition would violate 
unitarity, as can be checked easily by taking the time derivative of the norm of the field. 

From a physical point of view the problem lies in the fact that the MIT bag model sets 
the field to zero outside the bag already in the Lagrangian. In the Bogolioubov model with 
a finite square well potential, there is a non-zero field also outside the well. If the wall 
moves inward, the field gains enough energy so that parts of it can go out of the well. The 
higher the potential the more energy is transferred to the field during a compression, mainly 
because the Dirac field at the wall does not approach zero as the potential goes to infinity. 

The MIT bag model overcomes this problem by means of the boundary condition Eq. [1^, 
which sets the field to zero at the moving boundary. The drawback, as we have seen, is the 
absence of a bag-like solution for a dynamical boundary. 

We also asked ourselves whether giving a mass to the fermions would alleviate the prob- 
lem. The static bag wavefunctions, as their massless counterparts, do not equal zero at 
the boundary. A moving boundary, however, requires also in the massive case the bound- 
ary conditions Eq. |l^ and hence, as mentioned before, the solution cannot be written as a 
time dependent combination of static cavity eigenfunctions. We could not find the general 
solution in this case, and we conjecture that no solution exists for massive fermions. Our 
conjecture is based on the fact that Eq. |12| implies a delicate cancellation of all the Fourier 
components of the wavefunction. In order to maintain this zero boundary condition at all 
times, all the Fourier components of the wave should travel at the same speed as that of the 
moving wall. This is indeed possible in the massless case, if the bag wall expands with the 
speed of light. For massive fermions however, each Fourier component travels at a different 
speed, and it may not be possible to satisfy the boundary condition Eq. [l^ at all times. 
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In this paper we have considered specifically the MIT bag model, which is singular in the 
sense that no kinetic energy term is associated with the motion of the boundary. However it 
is straightforward to see that our result applies without modification to a non-singular model 
such as the "Budapest" bag [0] which includes a surface tension besides the volume energy 
term. This is because the presence of a boundary kinetic term does not modify the linear 
boundary condition Eq. |, which is responsible for the unphysical solution. Analogously to 
the MIT bag, energy conservation would require the "Budapest" bag to have a zero surface 
tension, besides B = 0. 

In summary, we have shown that the equations of motion derived from the dynamical 
MIT bag Lagrangian with massless quarks and a spherical boundary admit only one solution 
corresponding to a bag expanding at the speed of light. This result raises the question 
of whether the quantization of the theory can provide stable solutions. We infer that a 
consistent dynamical bag model for absolutely confined fermions must include an interaction 
term with some other fields at least at the boundary of the domain. For example, in a work 



to be published elsewhere |]TT|] we will show that a dynamical chiral bag model |Q admits 
physically meaningful solutions. 

We thank the support of a Hong Kong Research Grants Council grant CUHK 4189/97P 
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